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ABSTRACT 


The advent of the missile firing submarine has added yet another 
dimension to the problem of defending convoys and task groups during 
ocean transit. The specific situation wherein the submarine must 
surface to fire a relatively short range missile against a convoy of 
ships is considered. The model developed considers several different 
problem parameters. It enables the calculation of probability of 
detection of the submarine, probability of killing the submarine before 
a particular missile is fired, and the expected number of missiles that 
the submarine will fire. Selected results from randomly selected 


parameter values are also presented. 
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CHAPTER I 


INTRODUCTION 


1.1 Objective. 


The constantly evolving process of measure and countermeasure have, 
in recent years, added yet another weapon to the submarine arsenal. The 
ability of the submarine to surface and launch an anti-shipping missile 
allows the attacking submarine to stand off many miles from its intended 
target and launch its attack. This same type of attack is also very 
possible against either strategic or tactical land based targets that 
are close to shore lines. 

In order to assist operational forces in gaining an insight into 
the problem, and to make guidelines available that will assist in the 
selection of numbers of aircraft, altitudes to fly, etc., a mathematical 
model is presented in this paper. This same model can be used with equal 
facility by prospective developers of future search systems. The =e 
presented is an adaptation, through set and probability theory, of 
search theory techniques developed by P. A. Morse, G. E. Kimball, and 
B. 0. Koopman. [2. | 

Included in this paper, besides the model, are probabilities of 
detection and kill, and expected numbers of missiles fired for randomly 


selected values of the various necessary parameters. 
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CHAPTER II 


MODEL DEVELOPMENT 


2.1 Assumptions. 


Several assumptions are basic to the development of this model. 

Some are made in light of the author's operational experience, and others 
are assumptions of necessity in order to release the model from human 
factor encumberances that cannot be readily expressed mathematically. A 
slightly different developmental path for the same model is presented in 
Appendix I, including the assumptions made and its mathematics. 

The basic scenario is described as follows. Consider a convoy 
(task group) at sea. It is known that enemy submarines present a threat, 
but the presence of one or more of these submarines in the immediate area 
is not an established fact. Among the offensive threats that the enemy 
submarine force has, is the ability to surface and launch anti-shipping 
missiles. An airborne radar search is to be conducted about the convoy 
as a defensive measure. 

The launching of a missile requires the submarine to remain on the 
surface for some finite length of time. The position of the submarine 
| within the area described about the convoy by the missile flight radius 
is unknown. It is assumed that this position is equally likely to be 
at any point in the laaeckiae area. Factors affecting the missile 
launch radius, hence the area to be searched, are discussed later in the 
paper, 


Two further assumptions about the submarine are made. First, it is 
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assumed that the submarine has no information as to the whereabouts of 
the search aircraft at the instant that he decides to surface. Second, 
the number of missiles to be fired is decided upon before surfacing, 
thereby fixing the length of the firing cycle before the actual time of 
surfacing. This latter restriction is not extremely unreasonable in 
the case where the submarine has orders bial shook at a specific time 
(i. e., the start of a surprise attack leading to open hostilities). 

The only search weapon considered in this model is the aircraft 
search radar. Further, it is assumed that the paths flown by the 
aircraft are random within the search area. A study of Figure (1) will 
show this to be a reasonable assumption, even when each aircraft is as- 
signed a search sector within the search area. Note that search sector 
lines could be drawn in the figure, but without them drawn the randomness 
becomes quite clear, 

A parameter called the "Assured Radar Detection Range" is used 
throughout —_ model, This quantity is defined as the range at which 
the number of targets tebea at lesser ranges equals the number of 
targets contacted at greater ranges. This has been called the "Cookie 
Cutter" detection range. From this definition, the probability of de- 
tecting a target passing within this range of an aircraft is ome, Care 
im selection of the "Assured Radar Detection Range" for a particular 
situation is to be ey This subject will be discussed later in 
the paper. 

Once the submarine has been detected it is assumed that the only 


possible attack vehicle is the detecting aircraft. If any other aircraft 
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Seale: 1 inch = 50 miles 
Search Parameters: 20 minutes at 150 knots, with 30 mile radar. 


Random Search Assignment Diagram 


Fagures 1 


™~ 
e8f tm be - ‘wont £ . 
~tabat stim OF diiw ,avond Oak te aeteoio OS Be ha "18" 


= es 


ii ? a pag wis 


meataeid JasangiseA dogse8 mobeeh || 


I eiupet 


were closer, this would imply that that aircraft would have made de- 
tection first. The latter assumption must be modified to some extent 
if more than one type of aircraft is present. In addition, the as- 
sumption is made that the firing of one or more missiles does not give 
away the position of the submarine, even though the presence of the 
submarine somewhere in the area is then obviously peeea: 

The delay involved in speed build up during the initial part of the 
attack phase, when the aircraft is transitioning from cruise to dash 
(attack) speed, is assumed offset by the ability of the aircraft to 
attack using short range standoff weapons such as air to surface rockets. 
Also, arriving in attack position implies a successful kill with proba- 
bility one for the purposes of this model. 

All other parameters used in this model are fairly straightforward. 
Either they are standard aircraft parameters, or they will have to be 
estimated, as ia the case of the enemy sitiea tide parameters. 

2.2 Probability of detection. 

The total area searched by each aircraft while the submarine is on 
the surface is depicted in Figure (2) or Figure (3), depending ee Eke 
sweep angle of the search radar in use. The calculation of the area 
involved is developed using: 

a. The area searched out immediately after the submarine surfaces 
is approximately equal to'the area covered by one sweep of the radar. This 
is some portion of the area of a circle, and is computed as Onin Where 


6 is the fraction of a 360 degree circle covered by the search radar, 
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and R is the assured radar detection range. 

b. In addition to the above the area covered while the submarine 
is on the surface is: 

(time on surface) X (aircraft cruise speed) X (twice the 
search radius measured perpendicular to the flight path) 

Cy gut 6 = 0.5 the search width (radius measured perpendicular 
to the flight path) is R sin COT» noting that 6b T=0 , iThis 
result can easily be seen by a study of the geometry of Figure (2). 

BE 6 = 0.5 then the search radius, as shown in Figure (3), is simply 
R. 

d. The time spent on the surface by the submarine will be de- 
pendent upon the time it takes to fire each missile, and the number of 
missiles to be fired. This time is computed using the formula 
(Tl + (m - 1)T2). Where Tl is the time from surfacing to launching of 
the first afesitet and T2 is the interval between successive launches. 
The quantity n is the total number of missiles to be fired. 

e. Combining the above subparagraphs together the total area 
searched by bes aircraft is then: 

0 =< O<=0.5 
= ON + al + G - 1)T2)V_R2sin( OTT ) (1) 
6.5 =O 100 
X= OTR? + (rl + (n - 1)T2)V_2R (2) 
To obtain the probability of detection of a submarine by N search 


aircraft the development used by P. A. Morse and G. E. Kimball is 
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followed. [2] In this development the random search postulated previ- 
ously is essential. In addition any overlap of search effort is assumed 
to occur in a random manner. 

Consider the instant when the total area searched since the subma- 
rine is surfaced is 2. The probability that the submarine is not detected 
prior to this instant is denoted by po, B). Next, increase the area 
searched by af. This increase being positioned at random within A, 
the total area to be searched. The probability that the submarine will 
not be detected in af is l - af /A. The probability that no contact 
is made while B + af is being searched is then denoted by po, {3 = dQ). 


Then: 


Po, B+aQ) “ [. - $f) po, 3) 
po, +d - P(O ieee 
A 


d P(O me) 


P(O, ) A 
Yielding the differential equation: 


aro, 2) =- #8 v0,Q ) (3) 


Solving this differential equation, with the boundry condition P(0,0) = 1, 


= Bia 


po, 2) ) =e (4) 


Therefore, since the total area searched in the actual problem is denoted 


results in: 


nO, the probability of detection of a submarine during its total on 


surface period is: 


_ NO 


PO) =10<6 °* (5) 
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2.3 Probability of detection of at least one submarine in a group. 


There are three possible surfacing situations for the multiple 
submarine problem. First, consider the case where two (or more) 
submarines surface simultaneously, As will be further elaborated on 
later, this would seem to be the most likely of the three cases. The 
probability of detection of at least one of M submarines is: 


P(detection of at least one of M) = 


i n°, _ nO ,™ 
= 2 

= 1.0 - 10-e e | 

0 

7 AO 

| A 
=1.0- le 

_ ™mOC 

A 

#=1.0- e (6) 


The second case to be considered is the situation where the subma- 
rines surface at disjoint intervals (no two or more on the surface at 
the same time). This is clearly the original problem stated and may be 
handled with the theory initially developed. 

The last case involves different surfacing times and overlapping 
surfaced intervals. Probabilities for this situation are calculable 
only for each individual case. In addition, this situation does not 
seem as operationally appealing as the previous cases. In this case 
the detection of one submarine would only serve to alert the convoy as 


to the presence of submarines in the area, hence making succeeding 
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attacks by succeeding submarines more risky to those submarines. 
2.4 Probability of kill before the X'" missile is fired. 

Now that the probability of detecting a surfaced submarine is calcu- 
able, the next logical question concerns what can be done about the de- 
tected submarine. Especially, what can be done before the submarine can 
fire his missile(s). 

In order to develop this part of the model a digression is necessary 
for Rackirtunld purposes. This development uses the methods demonst rated 
by E. Parzen. [3| 


Consider four events A, B, C, and D. 


P(A(\(BUC){)D) = P Ff ansyU aN] Md} 
=P [canBND)U (aNeNp)| 
= P(ANBMD) + P(ANCND) - PCANBNMDNal tcp) 
P(AMBMD) + P(ANCND) - P(ANBNCNMD) ~——s (7 
Next let: 


A = The event the submarine is in some searched area gg desig- 
mated detection in the following. 

B = The event that the submarine is in the initial sweep of some 
aircraft's radar upon surfacing, designated () in the following. 

C = The event that the submarine is fn the balance of a particular 
area of size io. designated [*| in the following. 

D = The event the submarine is killed, designated kill in the 
following. | 


Using the above definitions the following results are to be noted. 
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a.” BUC = ©+|-] vs) 
b. Bf)c s D (empty set) 
Using result (b), term number three of equation (7) can be handled thusly: 
P(AMBIICMD) = KAMDMD) =P) = 0 (8) 
Using the above definite fone and the result shown in equation (8), 
equation (7) can be rewritten in the form: 
P [ detection () (detection in © U detection in [-] yN ki11| 
= P(detection[) detection in ©) M ill) + 
P(detection () detection in [-] a) kill) 
= P(detection)P(detection in 0) | detection)¢ P(kill | detection [) 
detection in ()) + P(detection)+P(detection in [-] | detection)-« 
P(kill | detection {) detection in [-]) (9) 
The following results concerning equation (9) follow directly from the 
previous definitions. 
a. detection {) detection in @ = detection in © 
b. detection f) detection in [] = detection in [:] 
c. detection f/\(detection in © U detection in [] ) = detection 
Using these results, equation (9) can then be rewritten in the form: 
P(detection {) kill) = P(detection)> Z (10) 
Where: 
Z = P(detection in © | detection)-¢ P(kill | detection in (@ ) 
+ P(detection in [-] Wetecttion)aP(ik ld | detection in [-] ) (11) 
Since: 
| P(kill [) detection) = P(kill | detection)-P(detection) (12) 


from the laws of probability, equation (11) or Z is actually the 
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P(kill | aeccceton): 

In order to uncondition equation (11) we again turn to probability 
theory. Consider the two events F and G. 

P(F) = P(F/)G) + P(FN Gc‘) 

: = p(F | G):P(G) + P(F | G°)-P(G’) (13) 
Where G° is defined as G complement. Next we make the following defi- 
nitions. 

F = The event the submarine is killed, designated kill in the 
following. 

G = The event the submarine is in a searched area, designated 
detection in the following. 
We then arrive at: 

P(kill) = P(detection)+P(kill | detection) + P(no detection)- 

P(kill | no detection) (14) 

From a practical'standpoint the P(kill | no detection) must be equal to 
zero. Thesetece: 

P(X) = P(kill before the x missile is fired) 


= P(detection before the x 


missile is fired)«Z (15) 
Where Z is as defined in equation (11). 

The next step in the development is to explicitly define all of the 
terms of equation (15) in terms of the measureable parameters. The proba- 
bility of detection term is similar to that expressed in equation (5), 


with the only change being the substitution of the yee missile for n in 


equation (1) or equation (2), whichever is appropriate. 
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Examining equation (11) we find the first term of Z to be 
P(detection in 6) | detection). This is equivalent to the event that the 
surfacing submarine is within range R of an aircraft when it (the subma- 
rine) reaches the surface, given the fact that the submarine is in one 
of the searched areas B. (See Figure (2) or Figure (3), whichever is 
appropriate.) Since the submarine is equally likely to be in any part 


of [3 (see initial assumptions), this term can be represented as: 


area in 0) 


P(detection in (0) | detection) a total area in 


OTT Rr? 0<(O<0.5 
OTT R? + aev_sin(O7T) (TI +(X - 1)T2) 
= (16) 


OTT R? 0.5=O=1.0 
OTT R? + 2RV_(T1 +(X - 1)T2) 


The next term to be considered is P(detection in Elig@erectveaye 
This is clearly 1.0 - (the applicable portion of equation (16)). 

The third term to be considered is the P(kill | detection in © ). 
The given event in this probability statement implies that the submarine 


is detected as it surfaces, and will have a period of time equal to 


(Tl +(X - 1)T2) remaining on the surface before the oh missile is fired. 


The distance that the aircraft can travel in this period is then 

(71. + Ge - 1)T2) Var If this distance is greater than or equal to R the 
aircraft will always be able to get to attack position. Therefore, 
P(kill | detection in @) will be one in this situation. 


If the distance the aircraft can travel is less than R there will 
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be some percentage of contacts that, even though detected, will not come 
under attack. Again using the assumption of uniform distribution of 
contacts, the circular uniform distribution applies within the circle, 
and is of the form: 


debt 2071 rdr _ (2r)dr 
r ATT r2 2 


R (17) 


(Tl + (x - 1)T2)v, 
F (TL + (X - 1)T2)V,) = (2r) dr 
2 
: R 


é +(x - ma) : 


R (17a) 


Combining the information from the preceding two paragraphs: 


(T1 + (X = 1)T2)v,)"  (T1+(X-1)T2)v, = R 
ies 
P(kill | detection in @)) = (18) 
1 (1+ & - 1)T2)V, =R 
The last term needed is P(kill | detection in [] ). Detection in the 
square itself implies two restrictions. First, detection will be made at 
range R, and second, detection is made after some period of time has 
elapsed since the submarine surfaced. If (Tl + (X - 1)T2)V, is less than 
R the aircraft will never be able to attack the submarine before missile X 


is fired. Therefore, the P(kill | detection in [J is zero for this case. 


LE (TL + KK = 1)T2)V, is greater than or equal to R, then there will 
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be some period after the time of initial surfacing that the aircraft will 
be able to get into attack position. Again under the assumption of 
equally likely submarine position in the square, there exists equally 
likely detection times given the fact that @etecetan wil be made. 
Defining T, as the elapsed time the submarine has been on the surface 


when detected (a uniform random variable), we have: 


' = Ee se = a = 
me ) (Tl + (XK ~ 1)T2) O=t' = (Tl + (X 1)T2) (19) 


Examination of the physical situation indicates that an attack will be 
possible only when the total on surface time minus the elapsed time on 
the surface when detected is greater than the time needed to get to the 


target for attack. This is formulated as: 


[or x». 172) T, | v,2R (20) 

(Tl + (x - LT, - RET,V, (20a) 
Vv, + (X - 1)T2) - R 

0< T= I (20b) 


The period between the upper and lower bounds on Ty. expressed in re- 
lation (20b) could be called the useful time for detection. Using 


relations (19) and (20b): 


V,(Tl + (XK - 1)TQ- R V,(T1 + (X - 1)TQ-R 
P re S$ \-", | 
d t d 
v,@71 +o(xe= T)T2)'t- "R 
Na 
m conbiperaciole sp aie virgtonINo Ss Sew, | 
(Tl + (X = 1)T2) d 
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Multiplying and combining terms, and defining: 


sin(OTT) op SOr= 015 


(22) 
1 0.5 =O <1.0 


Q= 


this development shows the probability of kill before missile X is 


fired to be: 


| - 86 71R° + 2RV (T1 + (xk - me 
l-e 


6 Wert + & - 1)12)?v, 


7) (23a) 
OTR + 2RV (TL + (x - 1)T2)Q 


- § Orr’ + 2RV (T1 + (x - ies 


2RV_.Q 


1 (23b) 


COMR + 2v (Tl + (X - 1)T2)Q)V, 
Equation (23a) is used when (Tl + (X - 1)T2)V, = R, and equation (23b) 
is used when (Tl + (X - 1)T2)V, >R. 

2.5 Expected number of missiles fired. 
In order to be able to make this calculation the probabilities of 
the submarine firing exactly one missile, exactly two missiles, 

» through exactly n missiles are needed. At the present stage 
the model has yielded the probability of killing the submarine before 
missile X is fired. The event that exact ly one missile is fired is 
disjoint from any other firing event, as is the event that exactly two 
missiles are fired disjoint from any other firing event, etc. With 


these considerations in mind we proceed as follows: 
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P, (0) = P(O missiles are fired) 
P, (1) = P(O missiles are fired, or one missile is fired) 


(24) 


Pn) = P(O or l or 2 0r.... .. or n-l missiles are fired) 
Then using the set of equation (24): 

PY (2) - PY (1) = P(one missile is fired) 

PY (3) - PY (2) = P(two missiles are fired) 


° (25) 


PY (n) - P,(n-1) = P(n - 1 missiles are fired) 
The set of equations (25) can be written in the general form: 
P, (i) = Pedi +1) - PY (i) Lal, tues oi Pe: ce | (25a) 
The probability of the submarine firing exactly n missiles can then 
be computed by summing the (n - 1) equation (25a)'s and subtracting the 
result from one. 
n-1 


P,(n) =l - » PY (i) =l- PL (n) (26) 


The expected number of missiles fired may now be formulated in the 


normal way for discrete probabilities. 


E(x) = 0 P, (0) +1 P,(1) +2 P,(2) = pet Ra AS oi P,(n) 
= 0 PY (1) +1 (P,(2) - P,(1)) oc. ea dW = P(n)) 
=n- P, (1) - P.(2) a een eae ee P,(n) 
n 
=n- P_ (i) (27) 
K 
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CHAPTER III 


CONCLUSIONS 


3.1 General results. 

Representative values of probabilities of detection, probabitities 
of kill, and expected numbers of missiles fired are presented in 
Appendix B. The parameters used in this appendix were chosen at random 
from the many combinations of numbers used during development of the 
model. All computations were done using a digital computer. 

The manner in which the probability tables are presented in the 
appendix clearly shows two very interesting results of the model. The 
first of these is evident when, for a given set of fixed heading 
parameters (aircraft and submarine capabilities, and area to be searched), 
any range R is picked. Following across the line that was chosen in 
the tables it is evident that the return per added aircraft, in terms 
of increased probabilities, diminishes as the number of aircraft is 
increased. This of course is a characteristic of exponential curves, 
and is to be expected here in view of the development. 

A study of the PY (X)'s, for X equal one and two especially, illus- 
trates the second of the above results. nein using fixed heading 
parameters, but this time holding N constant and changing R, this study 
reveals that in some cases the probabilities of kill start to increase 
with increasing R, reach a maximum, and from there on will aindetett 


In order to explain this phenomenon, consider the circular shaped 
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area that the search aircraft's radar sweeps out during each sweep. When 
R is increased, the area searched is increased as a function of R squared, 
hence significantly increasing the total area searched. This also in- 
creases the probability of detection. If the aircraft could reach, and 
attack, all of the targets it detects before the missile of interest can 
be fired, this would imply a steadily increasing P(X). This is in fact 
the case for reasonable, tactically interesting sets of parameters. On 
the other hand, if R is extended past the distance of possible travel by 
the aircraft before the missile of interest is launched, two events take 
place. First, the probability of reaching all targets detected in the 
initial sweep of the radar after the submarine surfaces is now less than 
one. Second, the aircraft will be unable to reach any targets detected 
after the surfacing time in order to stop the missile launch of interest. 
Thus, considering the rate of increase in searched area with an increase 
in R, there exists a higher and higher probability that a submarine de- 
tected will be outside attack range. The above facts then give an 
intuitive reason why the P(x) 's do reach maximum values. After testing 
many sets of parameters that could be tactically interesting, the author 
found that the maximum occurred at the point where R = (Tl + (X-1)T2)V_. 

A partial proof of the above will be sketched through a combination 
of numerical and analytical techniques as follows: | 

a. Proof that equations (23a) and (23b) are in fact equal at 
R=(Tl + (X - 1)T2)V_> thus establishing continuity at this point. 

b. Proof that equation (23a) is monotone decreasing for: 


(Tl + (X - 1)T2)V, = hi— Co; (28) 
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c. Statement of a result about the form of equation (23b), that 
was observed from many numerical results, for: 
O<R = (Tl + (X - 1)T2)V,- (29) 
At this point, in order to simplify the notation, the quantity T will be 
used to refer to (Tl + (X - 1)T2). 
A comparison of the two parts of equation (23) shows that step (a) 


above requires that: 
2.2 
OTTT'v, 2RV_Q 


= l - 
Or + 2RV_TQ CAITR + 2V, TOV, 


at the point where R = TV,. Therefore, upon substitution we get: 


d 
ATR eK OTTR + 2v 79) % - 2Rv 79 
OTR? + 2RV TQ ( O7TR + 2v 19) & 
OTR? O7TR? (30) 


OTT? + 2RV_TQ Oj + 2RV_TQ 


The proof of step (b) is effected by examination of the first 
partial derivative of equation (23a) with respect to R. After some 


algebraic manipulation this is: 


N 2 
P(X) Nc A 70 Rr? +2Rv_7Q) 
se wx [1+ O 70x? + 2nv,20 Ege Se S 
y : (31) 
Moers: aC 6 TU x2 + 2Rv 79) 
Cc 2 2 
e OMe + 2.10) 1 T¥,.) 
K = (31a) 


(union: 2RV_TQ)7 


This derivative must be less than zero for all values of R greater 
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than TV, in order for the slope of the P(x) curve to be negative in 
this range, Examination of the constant term K shows that this is 

always a positive constant. Therefore, the term in the brackets in 
equation (31) must be negative for this proof to work. A series ex- 


pansion of the exponential term yields: 


N 2 
Nc ATT R* + 2Rv_7Q) 
ae ¢ -1+5¢6 71 Rr” + 2Rv 79) + 


BCO 70 kr” + 2RV TQ) 7 
yp A ee ces 


= sb pet) 


Substituting equation (32) into equation (31) and performing the sub- 


traction leaves: 


bP) (Rc O71 R? + 2ev 799]? Ars OTT Rr? + 2Rv 19)}° F 
OR nibh sa sa | 
+ exert: (33) 


Equation (33) is clearly negative for all allowable values of R under 
consideration, thus completing the proof of this step. 

Lastly step (c) is somewhat more difficult. Again a study of the 
first partial derivative of equation (23b) would seem to be the desired 
path of approach. The idea being to prove this derivative positive for 
the desired range of R. This derivative is (again after some manipulation): 

BF) 

yR rs 


xX Ou wo 2RV TQ) 
COTTR+ 2V.TQ)V, e 
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(34) 
6 TTR + av ta 


Equation (34) is quite complicated, and any definitive analytical re- 


l-e 


sults are very hard to obtain. The equation will break down in many 
forms, some parts of which are clearly positive, other parts clearly 
negative, but in no case found were any of the terms found sufficiently 
alike so that a statement about their relative magnitudes could be made. 
In view of this difficulty many computer calculations were made. Many 
widely varying combinations of parancters were used in these calcu- 
lations. All computations showed that equation (23b) did reach a 
maximum, and then start decreasing. In all cases where tactically 
interesting sets of parameters were used this maximum was found to be 
at ranges greater than TV so that when the maximum was reached this 
equation did not apply to the problem. 

The considerations involved in selecting "tactically interesting” 
parameter combinations were avoidance of unreasonably long on surface 
times for the submarine, use of so many aircraft and such large R's 
that NO( was many times larger than A, etc. 

The desired conclusion then is that by using reasonable sets of 
parameters the maximum P(X) 's will be at ranges equal to TV,- Tf 


unusual aircraft parameters are to be used care must be exercised in the 


selection of areas to be searched, etc., in order to derive meaningful 
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results from the model. 
3.2 Parameter relationships. 

Following up on the discussion in the latter part of the previous 
section, one of the most striking parameter relationships occurs in the 
PY (x) equations. Of particular note is the fact, that when the killing 
of the submarine by the detecting aircraft is desired, there is no bene- 
fit in increasing R above the range equal to (Tl + (X - 1)T2)(V,). In 
fact, ground is lost when R is greater than this limit. This has an 
affect on search altitude selections, and radar scope range settings, 
since these factors directly affect R. Furthermore, the designer of 
future aircraft has a handle on the relationship between R and Va for 
this type of problem. 

The area to be searched, when circular in shape, has a tendency to 
become very large, even with seemingly small increases in circle radius. 
Again, the area is a function of the radius squared. An unnecessary 
increase in the area to be searched can drastically reduce the effectiveness 
of the total search. 

‘Noting that the increase in the area actually searched by each 
aircraft is a function of both R and R squared, it becomes apparent that 
for a plain search problem the more range the better. An increase of 
just five miles, say from 20 to 25 miles, makes a 56% increase in the 
area covered by one sweep of the radar. 

3.3 Tactical implications. 
While this model enables the calculation of probabilities of 


detection and kill, and the expected numbers of missiles fired by the 
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submarine for specific sets of parameters, this is clearly not the whole 
picture. 

As an example, the selection of the size of the area to be searched 
is quite important. Should the radius of the circle in question be equal 
to the maximum missile travel distance? This author submits that the answer 
to this question is no. Besides this maximum launch range the planner 
should take into account possible range maximums at which the submarine 
might detect the convoy, disposition of the ships in the convoy, estimates 
of the reliability of the enemy missile (if available), possible other 
methods of attacking a detected submarine, among other factors, By 
making rational estimates of the above factors the search probability can 
probably be raised, as might be the kill probabilities. 

The “Assured Radar Detection Range" is a key quantity in this model, 
hence a good understanding of how values for it are arrived at is neces- 
sary, Many variables are involved. These include sea state, operator 
training and proficiency, proximity of other small non-submarine targets, 
flight altitude, meteorological conditions, and type of radar installed 
in the aircraft. It is stressed that this quantity is the range at 
which the total number of targets missed at lesser ranges equals the 
total number of targets detected at greater ranges. Care should be 
exercised in order to avoid overestimating this parameter. 

3.4 Model extensions to similar situations. 

Two situations that bear directly on this problem were not con- 

sidered in this development due to their necessary introduction of some 


human factors variables. The first is the case where the submarine delays 
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surfacing due to his ECM detecting a search aircraft's radar. The other 
is the case where the submarine surfaces, fires part of his desired 
launch of missiles, and then decides to submerge for the same reason 
enumerated immediately above. Both of these cases are very real possi- 
bilities, and consideration should be given them in the planning phase 
of a search operation. Since the ECM detection range of a radar is 
greater than the detection range on the submarine by the aircraft radar, 
these cases have some interesting ramifications. Among these is the 
calculation of the area covered by electronic emmissions that are of 
such strength as to dissuade the submarine from surfacing for some 
given period of time. This type of situation could delay missile 
launches, hence causing problems in the fire control solution. Also, 
using a random search will cause the submarine to be more cautious, 

even after the present aircraft passes (assuming an ECM detection), 
since the submarine does not know when that aircraft may turn and 
recover some given position, or even another aircraft might pass 
through the same area. Detailed extensions in these areas would be of 
great value, but are beyond the scope of this paper. 

Another obvious extension of this model is to the case of a search 
for a missile launching submarine on station off a aouet line. Justi- 
fication would be very difficult for the assumption of a uniform distri- 
bution of targets along an entire coast. This problem can be resolved 
by selecting areas with high probabilities of containing submarine 
missile launch stations, and then assuming a uniform distribution of 


targets throughout each individual area. 
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Applications of the basic probability of detection model presented 
can be extended to the search for trawlers, Search and Rescue operations, 
or for any situation that meets the assumptions made in the development 
of the model. When considering targets that are continuously on the 
surface, then becomes the total area searched by each aircraft while 
on station. 

Last to be considered is the case of multiple types of aircraft 
involved in the same search. In this case an({ must be calculated for 
each aircraft, based on its parameters. These (X's are then summed. 
The probability of detection computations are then straight forward. 
Some problems are encountered in the probability of kill computations. 
There exists the chance that an aircraft with a long range radar might 
detect the target, and there is an aircraft with a shorter range radar 
closer to the target than the detecting aircraft, In this situation 
the possibility of mutual assistance must be considered. The various 
factors involved in equations (23a) and (23b) will also have to be 
weighted for the different types of aircraft. This could be done by 
using the circle and square technique of section 2.4. This latter path 
is necessitated by the possible differing dash speeds of the aircraft 


involved. 
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APPENDIX A 


COMPARISON OF TWO RANDOM SEARCH MODELS 


Another approach to the search phase of this problem is presented 
by B. 0. Koopman. [2] This development will be outlined in this appendix, 
and some interesting comparisons pointed out. 

Three assumptions are made in this development. They are: 

a. The target's position is uniformly distributed in the area to 
be searched. 

b. The observer's path(s) is (are) random in A, in the sense that 
it (they) can be thought of as having its (their) different (not too 
near) portions placed independently of one another in A. 

c. On any portion of path which is small relative to the total 
length of path, but decidedly larger than the range of possible de- 
tection, the observer always detects the target within the lateral 
range R, and never beyond. 

All parenthetical expressions in the foregoing, except for "not too 
near", have been inserted by this author for the purpose of clarifying 
analogies that will be drawn later. 

The model development itself proceeds by dividing the total path 
length (call it L) into m equal segments. The probability of failing 
to detect the submarine, during its entire on surface period, is then 


the product of the detection failure probabilities for each segment. 


P(fail to detect the submarine) = (1.0 - Oya (Al) 
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P(detection) = 1.0 - (1.0 - Qye (A2) 
As m increases equation (A2) can be approximated by: 
m OC 
A 


P(detection) = 1.0- e (A3) 


for large m. It is immediately noted that equation (A3) is precisely the 
result obtained in the previous formulation of the problem. 

It would seem intuitively appealing to let the number of segments m, 
be equal to N, the number of search aircraft in the problem. The desire 
being to tie m, in some simple manner, to some parameter in the previous 
formulation. One also desires that the turns necessary for each indi- 
vidual aircraft be kept to a minimum so that the probability of de- 
tection on each leg is as close to (X/A as possible. 

A problem arises with assumption (c) when large (40 mile plus) 
radar detection ranges are used. For example, if the submarine is on 
the surface for 20 minutes, and the aircraft flies 50 miles during this 
time, it is clear that assumption (c) is not fulfilled. 

A comparison between the results of the exponential and equation 
(A2) formulations was made using a digital computer for the computations. 
A selected portion of these results is displayed as Table Al. 

A few words of explanation about the selection of the particular 
values presented in Table Al follow. Many combinations of parameters 
were run through the computer, from which the results presented in the 
table are representative, It was discovered that all parameter combinations 


yielded the same general trends, which can be noted from the entries in 
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TABLE 1A 
Constant Parameters: 
Oisrave ¥. = 150 knots Tl = T2 = 5.0 minutes n = 4 missiles 
All areas are circles with radii = 50 (50) 200 miles 
First row is exponential result. Second row is 1 - (1 - oS result, 


1 2 3 4 5 6 7 8 9 10 


-738 .778 = .813 
-769 .808 #.840 


week OC/A > 1.0 
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the table. Where no entry appears in the table, that particular situ- 
ation was considered to be tactically uninteresting. This event occurred 
when the total area searched by all aircraft was considerably larger than 
the total area to be searched. In all cases though, at least the values 
for N equal one are given. 

Three definite trends are evident from the tabled results. 

a. Given A and N fixed, the P(D)'s differ by greater and greater 
amounts as R is increased. 

b. If R and N are fixed, the difference in the P(D)'s decreases 
as A increases. 

c. Varying A and R will produce convergence, divergence, or some 
combination of the two for the range of N that is tactically interesting. 
It should be kept in mind that the differences referred to above are 
potwesni bs two different methods of formulating the problem. 

Next some of the causes for these trends are explored. In the 
process a metric will be developed for when to use which form of the 
P(D) equation. 

The reason behind trend (a) is related to assumption (c). Note 
that as R increases this assumption fits the actual situation less and 
less. In fact, rather than the high probability of a straight course 
that is desired for each aircraft, as R increases the probability that 
an aircraft will have to turn increases. A study of Figure (1) illus- 
trates this point. 

| Trend (b) can also be explained by the lessened probability of turn 


for each aircraft as A increases. When an aircraft turns there must be 
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some overlap (i. e., less total area searched), hence the probability 
of detection on each non straight leg is less than O(/A. 

The convergence/divergence problem in trend (c) is a function of 
the relative sizes of C(/A and N. The major point to note is that in 
no case, considering the range of N that is tactically feasible, is a 
clear cut case of close convergence established. 

From all of the foregoing two important points should be retained. 
First, each time assumption (c) fails to hold a lesser probability of 
detection is implied. Secondly, the exponential formulation always 
provides this lesser probability. Im the case of large A and rela-- 
tively small N (i.e., turns by each individual aircraft are less 
likely), the agreement between the two formulations is very close. 

All of the above points to the exponential formulation for proba- 
bility of detection celculeriens: It is to be noted that the two 
developments ended up in essentially the same place, but that the 
latter method deviates somewhat from its assumptions, thus making the 


method somewhat questionable. 
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APPENDIX B 


SAMPLE RESULTS 
Selected results using randomly picked parameter values are 


presented. The first eight pages contain the results of P(D) 


calculations, and the latter pages give P(X) and E(X) results. 
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APPENDIX C 


COMPUTATIONAL AID 


This appendix, through the use of the computational flow sheets 
on the following pages, is designed as an aid in the calculation of 
results that may be desired, but that are not presented in Appendix B. 
Needed exponential and sine function values will be found in Appendix D. 
Some additional symbols are introduced for computational ease. 
Their use in each case is self-evident, and none of the new symbols 


is used outside of this appendix. 
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APPENDIX E 


COMPUTER PROGRAM 


The computer program included in this appendix is designed to 
make all P(D), PL (xX), and E(X) calculations mentioned previously. 
The exponential formulation is used. 

The program is written in FORTRAN-63, and is divided into sub- 
routines for the various calculations. If only P(D) results are 
desired, Subroutines KILL and EXPNR can be removed, and vice versa. 
Comment cards have been used for ease in tracking through the logic. 


The program was checked out on a CDC-1604. 
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